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Signature of the Fragmentation of a Color Flux Tube 


Cheuk-Yin Wong^’Q 

^Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 

The production of quark-antiquark pairs along a color flux tube precedes the fragmentation of 
the tube. Because of local conservation laws, the production of a q-q pair will lead to correlations of 
adjacently produced mesons (mostly pions). Adjacently produced mesons however can be signalled 
by their rapidity difference Ay falling within the window of \Ay\<l/{dN/dy), on account of the 
space-time-rapidity ordering of produced mesons in a flux-tube fragmentation. Therefore, the local 
conservation laws of momentum, charge, and flavor will lead to a suppression of the angular correla¬ 
tion function dN/{dA(j> dAy) for two mesons with opposite charges or strangeness on the near side 
at {A(j), Ay)~0, but an enhanced correlation on the back-to-back, away side at A(f>~TT, within the 
window of | Aj/|<l/(ciA'^/dy). These properties can be used as signatures for the fragmentation of a 
color flux tube. The gross features of the signature of flux-tube fragmentation for two oppositely 
charged mesons are qualitatively consistent with the STAR and NA61/SHINE angular correlation 
data for two hadrons with opposite charges in the low-pT region in high-energy pp collisions. 

PACS numbers: 13.86.Hd,1366.Be 


I. INTRODUCTION 

In the production of particles in the central rapidity 
region in high-energy hadron-hadron collisions and e“'"e“ 
annihilations, the low-pp part of the spectra falls within 
the realm of soft nonperturbative QCD physics and is 
usually considered to arise from the fragmentation of a 
color flux tube (or its idealization as a QCD string) [T]- 
miHin]. The high-pT part in hadron-hadron collisions 
is considered to arise from a relativistic hard scattering 
of partons and subsequent parton showering [ini l20H34] . 
In the case of high-energy e+e“ annihilation, the annihi¬ 
lation leads to the production of high-energy quark and 
antiquark partons which are then subject to perturbative 
QCD parton showering processes to lead to the produc¬ 
tion of hadrons. 

Recently, it was found that the hadron pp spectra 
spanning over 14 decades of magnitude from about 0.5 
GeV/c to the highest pp at central rapidity in pp col¬ 
lisions at LHC energies can be adequately described by 
a Tsallis distribution with only three degrees of freedom 
in a form phenomenologically equivalent to the 
quasi-power law introduced by Hagedorn [39] and others 
m for relativistic hard scattering. The simplicity of the 
Pp spectrum suggests that a single mechanism dominates 
over the domain with pp>0.5 GeV/c at central rapidity 
in these high-energy collisions. As the high-pr region is 
known to arise from the relativistic hard-scattering pro¬ 
cess [ini iMsi], one is led to the suggestion that the 
hard-scattering process dominates over the domain with 
Pt>0.5 GeV/c in these high-energy pp collisions. Ad¬ 
ditional experimental evidences have been uncovered to 
support such a suggestion [55] . 

The dominance of the hard-scattering process does not 
imply the absence of the soft flux-tube fragmentation 
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process. It only stipulates that contributions from the 
hard-scattering process increase with increasing collision 
energies and the fraction of the contributions from the 
flux-tube fragmentation process becomes smaller in com¬ 
parison, as pointed out earlier in |28l [29]. As a con¬ 
sequence, there will be a transverse momentum bound¬ 
ary pph which separates the lower domain of flux-tube 
fragmentation from the higher domain of hard-scattering 
dominance. 

As the pp collision energy decreases, the role of the 
flux-tube fragmentation and hard scattering will be re¬ 
versed, with an increase in the fraction of contributions 
from the flux-tube fragmentation and a shift of the trans¬ 
verse momentum boundary pp^ to greater pp values. It 
is of interest to see how the two processes interplay and 
how the boundary function ppb{y/s~^) between the two 
processes depends on the collision energy -^snn- In ad¬ 
dition to being an intrinsic physical property of the pp 
collision process, the boundary function ppbir/sj^) sep¬ 
arating the two processes in pp collisions may have im¬ 
plications on the early evolution dynamics, the thermal- 
ization of the produced medium, the quenching of jets, 
and the formation of the quark-gluon plasma, in high- 
energy nucleus-nucleus collisions. It is therefore desirable 
to search for ways to discriminate the process of flux-tube 
fragmentation from the process of the hard scattering in 
pp collisions so that they can be separated out and the 
boundary ppb mapped out as a function of the collision 
energy. We need well-defined signatures for the flux-tube 
fragmentation and the hard-scattering processes. 

Two-hadron AejyAfj angular correlations have been 
previously first suggested and used by the STAR Col¬ 
laboration to separate phenomenologically the ‘soft’ 
and ‘hard’ components in nucleon-nucleon and nucleus- 
nucleus collisions [511155] . The signature for the hard 
component, represented by the production of two back- 
to-back jets (minijets) in the hard-scattering process, 
is well known [TOl Esm [30], as indicated for exam¬ 
ple by STAR AcfhAr] correlation data for two oppositely 






2 


charged hadrons with pt>0-5 GeV/c in pp collisions at 
^/st^=200 GeV shown in Fig. 1(b) [4llfe] . It consists 
of (i) a near-side cluster of particles at (A^, A? 7)~0 for 
one jet, and (ii) an away-side ridge at A(/)~7r along Arj 
for the other jet. 
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FIG. 1. (Color online) The angular correlation data 
Ap/y/p^cf as a function of the azimuthal angle difference 
(j>A~A(j) and pseudorapidity difference rjA—Arj of two oppo¬ 
sitely charged hadrons in pp collisions at yG]Gv=200 GeV 
(from Fig. 3 of |42] of the STAR Collaboration). Fig. 1(a) is 
for hadrons in the domain with pt<0.5 GeV/c. Fig. 1(b) is 
for hadrons in the domain with pt>0.5 GeV/c 

On the other hand, the ‘soft’ component, defined 
phenomenologically by the STAR Gollaboration to be 
the mechanism that dominates the production process 
in the domain with p7’<0.5 GeV/c for pp collisions at 
V®AfJV=200 GeV, was associated with the {Acj), Ary) cor¬ 
relation for two oppositely charged hadrons as shown in 
Fig. 1(a) [HMdSj . The small narrow peak at (A^, Ary) ^ 
0 is an experimental detector artifact. After removing 
this artifact, the correlation is suppressed at (A^, Ary)^0 
but enhanced at (A())^7r, Ary~0). The two-hadron corre¬ 
lation of the soft component was also taken phenomeno¬ 
logically as a one-dimensional Gaussian distribution at 
A?y~0 independent of Acj) [4TV145) . A rigorous proof of 
the microscopic connection between the two-hadron A(j)- 
Ary correlation and the signature of flux-tube fragmenta¬ 
tion is still lacking. The other suggested signature of a 
one-dimensional Gaussian at Ary~0 also needs to be reex¬ 
amined carefully. Starting from the constituent level, the 
microscopic connection of dN/{dA(j)){dAr]) in a flux-tube 
fragmentation will be the subject of the present investi¬ 
gation. 

The search for the signature of the fragmentation of 
the flux tube is also useful for a better understanding 
of the multiparticle dynamics of the fragmentation pro¬ 
cess. Experimental verification of the signature will pave 
the way for further investigations to find out how a lead¬ 
ing quark and an antiquark pair can produce a chain 
of hadrons that are ordered in space-time and rapidi¬ 
ties through the nonperturbative QCD processes. Of 
particular interest is the possibility of an event-by-event 
exclusive measurement in which all of the momenta of 
the produced particles in a reaction event are measured, 
matched, and correlated, as in a jigsaw puzzle, to see how 
the space-time-rapidity ordering may allow the event re¬ 


construction of the chain of hadrons at the moment of 
flux-tube fragmentation. 

To search for the signature, we envisage that in a flux- 
tube fragmentation in hadron-hadron collisions or e“'"e“ 
annihilations, a flux tube is initially formed between the 
leading quark and antiquark (or diquark in the case of 
a nucleon-nucleon collision) when they pull apart from 
each other at high energies. The vacuum is so polar¬ 
ized that ordered pairs of quarks and antiquarks are pro¬ 
duced inside the tube via the Schwinger pair-production 
mechanism or the QED2 inside-outside cascade mech¬ 
anism [3 0 HMHl HO]- The interaction of the 

produced quarks with antiquarks produced in adjacent 
vertices leads to the production of mesons and the frag¬ 
mentation of the flux tube. 

The production of a quark-antiquark pair needs to 
obey local conservation laws. These conservation laws 
will impose constraints and will lead to correlations of 
various quantities, for two longitudinally adjacent pro¬ 
duced mesons. The occurrence of two longitudinally ad¬ 
jacent produced mesons, on the other hand, is signalled 
by their proximity in rapidity, on account of the space- 
time-rapidity ordering of the produced mesons (see for 
example, Ghapter 6 of uni). Hence, the local conser¬ 
vation laws will lead to correlations of various quanti¬ 
ties as a function of the azimuthal angle difference Acj) 
and rapidity difference Ay (or approximately, the pseu- 
dorapdity difference Ary) for a pair of produced mesons 
in a flux-tube fragmentation. These angular correlations 
can be used as the signature of the flux tube fragmenta¬ 
tion in high-energy nucleon-nucleon collisions and e“'"-e“ 
annihilations. 

In this paper, we restrict our attention to the central 
rapidity region. In Section II, for a system of four pro¬ 
duced constituents we single out the relevant degrees of 
freedom to describe the flux-tube fragmentation process. 
In Section III, we discuss the momentum distribution of 
two nonadjacent mesons when the quark and antiquark 
constituents in the detected mesons are produced without 
correlations. In Section IV, we examine the transverse 
momentum distribution of adjacent mesons in which the 
quark of one meson and the antiquark of the other meson 
are produced at the same point, subject to local conser¬ 
vation of momentum. In Section V, we discuss the ra¬ 
pidity distribution of produced mesons in a flux tube. In 
Section VI, we examine the charge correlation of two ad¬ 
jacent and nonadjacent mesons for a system with two fla¬ 
vors. In Section VII, we consider the correlation of charge 
and strangeness for two adjacent and nonadjacent mesons 
for a system with three flavors. In Section VIII, we 
present the two-particle correlation function for measure¬ 
ments with different combinations of meson charges and 
strangeness quantum numbers. In Section IX, we present 
numerical examples of theoretical correlation functions of 
dN/dA<j)dAr] for different charge and strangeness config¬ 
urations to provide signatures for flux-tube fragmenta¬ 
tion. In Section X, we compare the theoretical signature 
of the correlation of two oppositely charged hadrons with 
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STAR and NA61/SHINE two-hadron correlation data for 
pp collisions [4TVl49j . In Section XI, we present our con¬ 
clusions and discussions. 


II. PRODUCTION OF TWO MESONS IN 
FLUX-TUBE FRAGMENTATION 

Our objective is to search for the signature of flux-tube 
fragmentation by studying the correlations between two 
mesons produced in a flux tube. Starting with the micro¬ 
scopic constituent momenta, it is necessary to enumerate 
the relevant degrees of freedom before we can examine 
the correlations of the observed composite mesons. 

We consider the production of two adjacent mesons in 
a flux tube at the moment of its fragmentation, as de¬ 
picted schematically in Fig. where the leading quark 
pulls apart in one direction while the other leading anti¬ 
quark pulls apart in the other direction at high energies. 
As the leading quark and antiquark pull apart, a flux 
tube is formed and the vacuum in the tube is polarized. 
The flux-tube fragmentation process is initiated by the 
production of many ordered q-q pairs along the flux tube 
by the Schwinger particle-production mechanism or the 
inside-outside cascade mechanism, as discussed in Refs. 
[I11I4H201. The production of these q-q pairs takes place 
locally at spatial points along the tube. 
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FIG. 2. Schematic depiction of the production of quarks 
and antiquarks in two adjacent mesons P 12 and P 34 as the 
leading quark is pulling to the left and the leading antiquark 
is pulling to the right, with the production of pi(antiquark), 
P 2 (quark), p 3 (antiquark) and p 4 (quark) along the flux tube. 
Subsequently pi interacts with p 2 to form meson Pi 2 , while 
P 3 interacts with p 4 to form meson P 34 . 

We focus our attention in a small section in the central 
rapidity region in the center-of-mass system and consider 
the pair p 2 (quark) and pa (antiquark) produced at the 
vertex V in Fig. We shall use the particle label also 
to represent the particle momentum. After p 2 and pa are 
produced, p 2 interacts with the adjacent antiparticle pi 
with a confining interaction, leading to the formation of 
the meson P 12 . Similarly, pa interacts with the adjacent 
P 4 to form the adjacent meson P34- In the leading order of 
the confining interaction, the constituents of meson P 12 
and the constituents of the adjacent meson Pa 4 do not 
interacl 0 because the linearly confining interactions are 


^ There are however higher-order residual interactions between 
mesons, as discussed for example by Peshkin and Bhanot ESI 
and in Refs. 15111521 . 


screened by their partners. The formation of the nearly 
noninteracting mesons P 12 and P 34 leads to the fragmen¬ 
tation of the flux tube into mesons along the longitudinal 
direction, populating the longitudinal momentum space 
in the form of a rapidity plateau. 

The pair-wise two-body interaction between pi and p 2 
as well as between p^ and P 4 along the flux tube makes 
it simple to study the many-body problem of quarks and 
antiquarks in two-body basis [52 |55j . We can use Dirac’s 
constraint dynamics [531459] to separate out the coordi¬ 
nates of a pair of particles into their center-of-mass and 
relative coordinates in a relativistically consistent man¬ 
ner. 

For the pair of interacting particles with 4-momenta 
Pi and P 2 and rest masses mi and m 2 , we construct the 
total momentum P 12 and the relative momentum pi 2 , 


where 


and 


P 12 = 

Pi +P2, 

(la) 

P2 • P 12 

Pi ■ P 12 _ 

(lb) 

p2 

^12 

Pi p2 P2, 

^12 

Pi ■ P 12 

P 12 +pI-pI 

(2a) 

p2 

^12 

^P?2 

P2 • P 12 

P 12 +pI-pI 

(2b) 

p2 

^12 

^P?2 


are the projections of the momenta pi and p 2 along the 
direction of the total momentum P 12 . The inverse trans¬ 
formation is 


Pi 


P2 


Pi ■ P12 p 
—“ 2 —P 12 +P 12 , 
W 2 

P 2 • P 12 

p2 

W 2 


P 12 — Pl2- 


(3a) 

(3b) 


Under the interaction 4 >(a;_Li 2 ) which depends only on 
the time-transverse relative coordinate a;_Li 2 between pi 
at a;i and p 2 at a; 2 , namely. 


{xi- X 2 ) ■ P 12 p 

a;_Li2 = (xi- X 2 ) - ^2 --^12: 

^12 


(4) 


with the property a;_Li 2 • P 12 = 0 , the eigenvalue equation 
for the meson bound state with wave function ( 2 ^- 112 ) 
in the relative coordinate is [551455] 

[&2(Pf2;"2?,m^) -bp ?2 - 4’(a;_Li2)]|V’(a;_Li2)) = 0, (5) 
where 


b‘^{Pf2,ml,ml) = 


P12 -2P'l2{ml+ml) + {m\-'mlY 


( 6 ) 


and 


p2 _ TiX2 
-<12 — -'' 212 ) 


Mi2 = 



(7) 


b'^ + mf + J + 1712 - 


( 8 ) 
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The solution of the bound state equation ([^ then leads 
to the meson wave function?/’(a;_Li2) and the meson mass 
Equations and Q indicate that as a result of 
the interaction $(x_li 2) that exists between pi and p2, 
the momentum elements are related by 

dpiS{pl - ml)dp2S{pl - ml) 

= dPi2S{Pl2 — Ml2)dpi2S{b'^—pl2 — ^ix±12))S{po,12)- ( 9 ) 


In the context of the fragmentation of the a flux tube, 
the produced antiquark pi and the quark p2 line up along 
the tube as in a one-dimensional string at the initial mo¬ 
ment of meson production, and one can use the string ap¬ 
proximation to separate out the longitudinal and trans¬ 
verse degrees of freedom by approximating the Pi2 in Eq. 

by Pzi2 and the confining potential coordinate x±i2 to 
be longitudinal. The longitudinal 12 and the transverse 
degrees of freedom can then be separated. The solution 
of the bound-state problem takes care of the Pzi2 degrees 
of freedom. Upon separating out and subsequently in¬ 
tegrating over the variables of To,12, Po,i2 and ^2,12, the 
remaining relevant degrees of freedom in the momentum 
element for particles pi and p2 in the flux-tube fragmen¬ 
tation are 

—p^—dPT 12 dpTl 2 = dyi 2 dPT 12 dpT 12 J (10) 
-c/12 

where yi2 is the rapidity of the meson P12, and the 
subscript T represents the transverse (two-dimensional) 
components. 

Similarly, the particles pa and p4 can be described by 
a composite meson with total momentum P34 


^34= P3 +P4, 


P4 • T 34 P 3 • P 34 

T34= —^2-T3-^2-T4, 

4^34 4-34 


(11a) 

(lib) 


with a relative momentum P34 containing a bound state 
with the meson mass P^^ = M'^^. The inverse transfor¬ 
mation is 


P 3 P34 T-\ , /I o \ 

P3 = ^2 -f34 +P34, (12a) 

-^34 

P 4 • T34 

Pi = 52 - -P 34 -P34- (12b) 


Thus, for two mesons P12 and P34 in the prob¬ 
lem of flux-tube fragmentation, the probabil¬ 
ity distribution dN is in general a function of 
dN{yi2, 2/34, Pti2 , Pt34 , Pti2 IPT 34 ) • 


III. THE TWO-PARTICLE TRANSVERSE 
DISTRIBUTION FOR NONADJACENTLY 
PRODUCED MESONS 

The process of q-q production for two adjacent mesons 
P12 and P34 depicted in Fig. [^occurs in the neighborhood 
of the vertex V. We envisage that similar processes of q-q 
production occur independently at many other vertices 
along the flux tube, leading to the production of mesons 
populating the longitudinal momentum space in the form 
of a rapidity plateau. 

We can consider two of the produced mesons which 
can be related to each other in two different ways. The 
mesons P12, and P34 can be in a (longitudinally) adjacent 
state as in Fig. in which the quark p2 of P12 and the 
antiquark p3 of P34 are produced at the same spatial 
point V and are correlated. The two mesons P12, and P34 
can be in a (longitudinally) nonadjacent state as in Fig.|^ 
in which all quarks and the antiquarks of the two mesons 
are independently produced at different vertices and are 
not correlated with each other. We shall often label the 
state of two mesons by the superscript X with X=A for 
the adjacent state and X=N for the nonadjacent state. 



FIG. 3 . Schematic depiction of the production of quarks and 
antiquarks pi, p2, ps, and p4 in two nonadjacent mesons P12 
and P34, in a flux-tube fragmentation. The quark p2 and the 
antiquark ps are not correlated because they are produced at 
different vertices. 

In the production of the four constituents Pi,P2,P3,Pi, 
we can separate out the longitudinal and transverse de¬ 
grees of freedom to write the distribution as 


The relevant degrees of freedom in the momentum ele¬ 
ment for particles P3 and p4 in the flux-tube fragmenta¬ 
tion are 

dy3idPT3idpT3i- ( 13 ) 


^ In practice, the quark and antiquark have spins, and one needs to 
use the two-body Dirac equation in place of Eq. in for the eigen¬ 
value equation involving a color-Coulomb plus a linear confining 
interaction. The numerical solution of the relativistic two-body 
Dirac equation then leads to the meson state with the meson 
mass Mi 2 and the proper meson spatial wave function, as de¬ 
scribed in Refs. | 6 QI I61| and references cited therein. A nonrela- 
tivistic model of a meson as a q-q bound state is given in Refs. 
[6l|63]. 


dN = dNydNr- ( 14 ) 

The probability for the two-meson state to be adjacent or 
nonadjacent depends mainly on the magnitude of their 
longitudinal rapidity differences jAi/l in comparison with 
the inverse of the rapidity density, l/{dN/dy), which we 
shall take up in Section V. We examine in this section the 
transverse momentum distribution of nonadjacent me¬ 
son, dN^, and in the next section the the transverse 
momentum distribution of adjacent mesons, dN^. 

The quarks and antiquarks of the constituents of the 
mesons are produced inside the flux tube, the wave func¬ 
tions need to obey the boundary condition appropriate 
for the tube. They acquire a transverse momentum dis¬ 
tribution governed by the geometry of the tube with a 
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standard deviation (t„ related to the radius of the tube 
-^tube by 

h 

^-■ (15) 

-^tube 

In addition, the quark and antiquark pair are produced 
in a Schwinger pair-production mechanism, from which 
they acquire a transverse momentum distribution with a 
standard deviation aq given by [131 IHl IHl [20] 


CTq = 



(16) 


where k is the strength of the linear color-electric held 
between the polarizing leading quark and antiquark. By 
folding the transverse momentum distribution from the 
hux tube with the transverse momentum distribution 
from the Schwinger mechanism, the transverse momen¬ 
tum distribution ciNxi of each constituent p^i can be 
represented by a Gaussian distribution as 


dNn 1 /- dg exp{-|^} 

dPxi (V^cr«)^ J {V^(7qy ^ 

where 


{pn - g)n 
2^2 J 

(17) 


a^=al + al (18) 

Therefore, when the four constituents pi,p 2 ,Pz,Pi are 
produced independently without correlations for two 
nonadjacent mesons, the transverse momentum distribu¬ 
tion dN:^ of these particles is 


dN^=. 


dprp-^dpx2 ^Px3 ^PxA 

(-v/27rtT)® 



Pti Pt2 Pt3 7*^4 1 

2 cr 2 /• 

(19) 


Upon averaging over E, we have 
1 1 .. . 1 


dN^ = 


4(-\/^cr)4 2 


l-dA^ 


Emax Emin JS 


dE 


X / dPT12dPT12PT34dPT34 SXp ^ ~ 


p2 I p2 

^T12 ^ ^T34 


which gives 


diVf 


1 


dA()) Stt^ 


(23) 


Consequently, for nonadjacent mesons in the absence of 
any correlation between the produced quarks and anti¬ 
quarks, the transverse momentum distribution dN^ is 
independent of the azimuthal angle difference A(j) of two 
nonadjacent mesons. 


IV. THE TWO-PARTICLE TRANSVERSE 
DISTRIBUTION FOR ADJACENT MESONS 


For two adjacent mesons as shown in Fig. the trans¬ 
verse distribution dV^ in Eq. (19) in the last section will 


need to be modihed to become 


dNx = dNqdNqidNqii —=- dpxidpx2dPT3dPTi, 

(V27rcr„)8 

_r (Pti- g') ^ + (Pt 2- g) ^ + (PT3 + <7) ^ + {Pt4 - g") ^ 

xexpl I, 

(24) 


where the momentum kicks g, q', and q" of the Schwinger 
mechanism have the distributions 

dNq^= Jl dq^ exp{-^}, g^=g,g',and g" (25) 
(V27rcr,)^ 2(7^ 


For the most likely case in the production of mesons with 
light quarks with mi = 7712 = 7713 = 7714 , the transformations 


from {pxi,PTj) to {PTij,PTij) can be easily 
In terms of the composite transverse mo¬ 
menta and their relative momenta, the momentum dis¬ 
tribution (19) becomes 


dN^ = 


(V^cr)' 
1 


-dPT12dPT34dpTl2dPT34 


X exp <j-^ ( -k 2p^i2 


T34 


+ 2Pt34 If ( 20 ) 


The relative transverse momenta Pxi 2 and Px 34 can be 
integrated out to yield 



We represent Pnj by {pnj , 4’ij ) and introduce the az¬ 
imuthal angle difference A(j) and sum E, 


A(f>— (j)l2 — 034, E — 012 + 034. (22) 


The quark p 2 and the antiquark p^ in Fig. [^are produced 
by the Schwinger mechanism at the point V. Because of 
the local conservation of momentum at the production 
point U, the transverse momenta Px 2 and Pxs are cor¬ 
related and the correlation is expressed in Eq. (24) as 
a shift of their momenta Px 2 and p^^ in Eq. (19) to 
Pt 2 — g and Pxs + g respectively. The constituents pi 
and p 4 themselves are products of Schwinger pair produc¬ 
tion and their transverse momenta Pxi and Px 4 become 
Pti ~ ^cid pj,^ — < 7 ", depending on the momentum 
kicks g' and g". 

After integrating over g' and g" in the distributions 
dNqi and dNqii as in Eq. (fl7|), the momentum distribu¬ 


tion dN^ in Eq. (24) becomes 


dNx = dNqdNxiq), 

dNxiq) = r— . r— ^dpTidpT2dPT3dPT4 

y\Z7r(TJ (^YZTTlT^j 

.. f Pti , (Pt2 “ Q? + (Pt3 + Q? , Pt4 1 

2^1 


(26a) 


.(26b) 
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Our task is to separate out the relevant degrees of free¬ 
dom in Eq. (26b). The most likely case is the produc¬ 
tion of mesons with light quarks, for which we can take 
mi=m 2 =m 3 =m 4 . By using the momentum transforma¬ 
tions of Eqs. ^ and (12), we convert Pt 2 ^ Pt?,^ and 
Pta, to Pti 2 ^ Pt 3 Ai Pti 2 , and Pt 34 - Equation (26b) be¬ 
comes 


clNriq) = - -7dPT12dPT34dpTl2dPT34 

X exp{A + B + C + D + E} . (27) 


where we collect the quadratic function of Pxi 2 ia A, 


-P 


T12 




,(28) 


A 2 r ^ I ^ 1 r-^T ’12 

we collect the quadratic function of Pt 34 in B, 

u— 2 r ^ I ^ 1 ,-Pt 34 , +-^734 + 29 ,, , 

B - -PT34 [^ + -PT34 ■ + ^2 ] (29) 

we collect the remaining quadratic function of Pti 2 in 

C, 


C=-P: 


1 


T 12 : 


1 1 
^ 2 a 2 + 2a?J 


- P 


T12 • 




(30) 


we collect the remaining quadratic function of Pt 34 in 

D. 

« = + (31) 

and the last term quadratic in in E, 

2 q^ 


E = - 


2al 


(32) 


By completing the squares of A and B in Eqs. (28) and 
(29), the integration over Pxi 2 and ^ 7^34 can be carried 
out analytically, yielding a function of Pti 2 , Pt 34 , and 
q. After these integrations, we obtain 


dNiq) = 


X exp 


X exp < — 


1 


{27rna^ + aiy 


-crlPTi2 


dpTi2dP 


T34 


al 


+ - 




T34 


T12 

8 cr ?2 


Pti 2 ■ q P 


T34 


2al 


8 cr ^2 


20-2 (j2 (j2 


Af.(33) 

^4/ 


2al 


where 

1 _ 1 1 
2 o -?2 20-2 + 20-2 ■ 

For the remaining degrees of freedom, we define ^12 as 
the azimuthal angle between q and Pti 2 ^ and ((>34 as the 
azimuthal angle between q and Pt 34 - We form the sum 
and the difference of the azimuthal angles 

^4>= 4’12 — (t> 34 i S = 012 + </i 34 ! (34a) 

0i2=i(S + A0), 034 = ^(S-A0). (34b) 


We then get the distribution dNxiq) defined in the in¬ 
tervals —27r<A0 < 27r and Emin(A 0 )<E < Sniax(A 0 ) 
where 

Emin(A0)= -27r -I- |A0|, (35a) 

Emax(A0)= -|A0| -I- 27r. (35b) 

The normalized distribution per pair of mesons P 12 and 



FIG. 4. (Color online) The transverse momentum distri¬ 
butions as a function of the difference A(j} of the azimuthal 
angles between two adjacent mesons, dN^/dA4> , and two 
nonadjacent mesons, dN^/dA(j) . 


P 34 is 

dNTjq) 

dAtp 


X exp 


X exp 


(A0) = 


(^max ^min)j 


dE 


PT12dPT12PT34dPT34 
2(271)2(^2+ ag)2 


^12 

2 


-CT^Pt12 

, 9 

2 2 

, eri2 

+erqPT34 

, 9 

2 

\ 

2(t2o-2 

2 

+ 2 

2cr‘^al 

2 

/ 


T12 

8 cr ?2 


PT12 ■ Q _ Pt34 

8 cr ^2 


2cr2 80 - 2 ^ 20-2 

where the scalar product Prij ■ Q is 

Prij ■ q = Pxij q cos(0ij). 


^^-^1.(36) 


erf, 


(37) 


Here, 0 i 2 and ^34 in the integrand can b e exp ressed as 
a function of A0 and E as given by Eq. (34b). We can 
calculate dNT{q)/dAcj) by integrating out E,Pti 2 , and 
Pt 34 in the above equation, and we obtain dNT{q)/dA(j) 
for different values of q. We need the values of and 
ag to get the numerical estimate. For a flux-tube radius 
of i?tube = 0.6 fm and k = 1 GeV/fm, we get cr„ = 0.33 
GeV/c, (Tg = 0.25 GeV/c, and cr=0.414 GeV/c. 

One finds that the distribution dNT{q)/dA(j) depends 
on the magnitude of the momentum |q|/crg whose distri¬ 


bution is given by Eq. (25). Upon integrating the dis¬ 


tribution of q, we obtain the distribution of A0 for two 
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adjacently produced mesons P 12 and P 34 , 


a constant given by 


^ dNriq) 


f exp|-d^|_ 

dA(/i J (v^crq )2 2 ct 2 '' dA(/) 


(38) 


The solid curve in Fig.|^is the dN^/dA(j) result obtained 
by a direct numerical integration over q. The correlation 
is suppressed on the near side at A(/)~0, and is enhanced 
on the back-to-back away side at A 0 ~ 7 r. Such a result 
is expected because the quark p 2 and the antiquark 
are produced in opposite azimuthal directions due to the 
local momentum conservation. The associated adjacently 
produced mesons P 12 and P 34 should also have back-to- 
back azimuthal correlations. 

Along with the correlation for an adjacent pair of 
mesons, there is also the azimuthal angular correlation 
dN^/dAcf) for a pair of nonadjacent mesons in a flux- 
tube fragmentation. The angular correlation dN^ / dAcj) 
corresponds to the case where the quark p 2 and the an- 
tiquark p^ ar e no t correlated by the momentum q. It is 
given by Eq. (231 with dN,^/dA(j)=l/8n‘^. It can also be 


alternatively obtained by setting q=0 in Eq. (36). 


The above results have been obtained by applying the 
local law of conservation of momentum for a system of 
equal quark masses. As the constituent quark mass of 
a strange quark is only slightly greater than the con¬ 
stituent quark masses of the up and down quarks [S3], 
the qualitative feature of the above results regarding a 
back-to-back correlation for adjacent meson production 
may also be applicable to cases involving the production 
of strange mesons. 


dN KTpj-o 
dy niT ’ 


(39) 


where k is the string tension and rriT is the transverse 
mass of the particle mT=\/'rrP + p'^ where m is the mass 
of a produced meson (mostly likely a pion). Thus the 
presence of a rapidity plateau is an indication of the ap¬ 
proximate validity of the occurrence of the the q-q ver¬ 
tices along a curve of constant proper time. We showed 
further in Exercise (7.2) of Ref. |20| that in that case, 
the produced particles are ordered in space-time-rapidity. 
The ordering of produced particles means that in the 
center-of-mass system, particles with a greater magni¬ 
tude of rapidity |y| are produced at a distance farther 
away from the point of collision and at a time later than 
and more separated from the time of collision, as shown 
in Eig. From such a picture of flux-tube fragmenta- 



V. RAPIDITY CORRELATION OF PRODUCED 
MESONS IN A FLUX-TUBE FRAGMENTATION 

The last section gave the azimuthal correlation of two 
adjacent and nonadjacent mesons in a flux-tube fragmen¬ 
tation. How do we identify adjacent and nonadjacent 
pairs of mesons? We shall use the space-time-rapidity 
ordering to correlate the rapidities of produced mesons 
with the spatial locations on the flux tube at the moment 
of fragmentation. 

Casher et al. |5] showed that in QED2 when a quark 
parton pulls away from an antiquark with the speed of 
light, the produced dipole density of produced q-q pairs 
is a Lorentz-invariant function and the lines of constant 
produced dipole density of produced q-q pairs are hyper¬ 
bolas with constant proper times. It is therefore reason¬ 
able to assume as a first approximation that the space- 
time distribution of the q-q production vertices should 
depend only on the vertex production proper time, Tp^o, 
relative to the point of the onset of separation of the 
leading quark and antiquark. 

Previously in Ref. [16] and in Exercise (7.1) of Ref. 
[ 2 ( 1 ] , we showed that if all the pair-production vertices of 
a fragmenting string fall on the curve of the proper time 
Tpro, the rapidity distribution of the produced mesons is 


FIG. 5. (Figure taken from Figure 6.4 of [20]). The space- 
time trajectories of the leading quark q moving to the positive 
direction (longitudinal direction) and the antiquark mov¬ 
ing to the negative direction, with the production of q-q 
pairs lying on the proper time curve of constant Tpro- A pro¬ 
duced q at one vertex interacts with the produced q in the 
adjacent vertex to form a meson. The trajectories of mesons 
are indicated by thick arrows. The rapidities of the produced 
mesons are ordered along the spatial longitudinal axis, and 
in time. 


tion, we envisage that produced mesons are ordered in 
space-time and in rapidity, with mesons that are pro¬ 
duced adjacent to each other closest in their magnitudes 
of their rapidity values m- Therefore, mesons are pro¬ 
duced adjacent to each other if their rapidity difference 
Ay falls within a rapidity width w, 


|At/| < w, 


(40) 


where w is the rapidity per produced meson, 
1 1 


dN/dy {3/2)dNch/dy 


(41) 


One notes from Eq. (39) that the width w is inversely 
proportional to the particle production time Tp^o- As 
the particle production time cannot be sharp and has 














a distribution (see Fig. 2 of m for an example of the 
scatter plot of the pair-production vertex proper times in 
the Lund model), the width w would also have a distri¬ 
bution. The condition relating Ay with adjacently pro¬ 
duced mesons should also have a diffuseness described 
by an additional parameter a. Therefore, within the 
nearest-neighbor rapidity width |A?/|;<r(; in the central 
rapidity region, the produced mesons are adjacent to each 
other and the rapidity correlation for an adjacent pair of 
mesons is 


dN^ 1 

dAy ~ 1 -I- ' 


(42) 


Upon approximating the rapidity y as the pseudora¬ 
pidity rj, the above correlation for adjacent pair becomes 


dN^ 1 

dAr] ~ 1 


(43) 


The complete correlation consists of both the transverse 
and longitudinal correlations. From the transverse az¬ 
imuthal correlation dN^/dAcf) for adjacent mesons as 
given in Fig. we have therefore the two-particle an¬ 
gular correlation for an adjacent pair of mesons as 


dN^ _ dN^ dN^ _ 1 dN^ 

dA(j) dAr] dAcj) dAy 1-|- 


On the other hand, outside this nearest-neighbor ra¬ 
pidity width with lAyl^w, the two mesons are nonadja- 
cent. The condition of being a nonadjacent pair, outside 
the nearest-neighbor rapidity width with jATyl^ry can be 
specified by 1/(1 -I- Consequently, the two- 

particle pseudorapidity correlation for a nonadjacent pair 
of mesons in a flux tube fragmentation is 


dN^ 1 

dAy ~1 -I- ’ 


(45) 


From the transverse azimuthal correlation dN^/dAcj) as 
given in Eq. (231, we have therefore the two-particle an¬ 
gular correlation for a nonadjacent pair of mesons as 


dN^ dN^ dN^ _ 1 _ 1 _ 

dA(f) dAy dAcf) dAy i _(-e(“-S tt^ 


In experimental measurements, one selects a pair of 
mesons whose charges and heavy-quark quantum num¬ 
bers can be the same or opposite. Different selections 
will result in different adjacent and nonadjacent frac¬ 
tions, and consequently different linear combinations of 
dN^/dAfpdAy and dN^/dAcj) dAy. They will lead to 
different shapes of the correlation functions. It is neces¬ 
sary to find out the adjacent and nonadjacent fractions 
for different charge and heavy-quark configurations. 


VI. CHARGE CORRELATION IN FLUX-TUBE 
FRAGMENTATION FOR QUARKS WITH TWO 
FLAVORS 


In a flux-tube fragmentation, the production of the 
quark-antiquark pair p 2 and at the local point V in 
Fig. [^must satisfy local charge and flavor conservations. 
This means that p 2 and ps must be quark and antiquark 
with opposite charges and flavor quantum numbers. Such 
a local conservation will induce correlations on charges 
and flavors of adjacent mesons. For nonadjacent mesons, 
there will be no such correlations. 


We can investigate all possible charge and flavor con¬ 
tents of two produced mesons if they are in the adjacent 
or nonadjacent state. We study the case of quarks with 
two flavors in this section, and the case with three flavors 
in the next section. 


We examine first the production of two adjacent 
mesons and construct in Table |T] all possible meson pro¬ 
duction configurations in which the quark p 2 and the 
antiquark p^ have different combinations of charges and 
flavor quantum numbers. Here, Qij is the charge of the 
composite meson with constituents pi and pj. 


There are altogether eight possible cases in Table [T] 
The number of cases A^"^(Qi 2 ,< 334 ) and the associated 
probabilities P'^(Qi 2 , Q 34 ) in different charge configura¬ 
tions (Qi 2 ,Q 34 ) for two adjacent mesons in a flux-tube 
fragmentation are given in TablejlT] We note that because 
P 2 andp 3 have equal and opposite charges, two adjacently 
produced mesons P 12 and P 34 cannot have charges of the 
same sign as shown in Table [III 


TABLE 1. Quark and antiquark configurations for two adja¬ 
cently produced mesons in a flux-tube fragmentation where 
P 2 and P 3 are constrained by local charge conservation and 
flavor balance. 


Pi 

P2 

Qi2 

P3 

Pi 

Q 34 

u 

u 

0 

u 

d 

-1 

d 

u 

1 

u 

d 

-1 

u 

u 

0 

u 

u 

0 

d 

u 

1 

u 

u 

0 

u 

d 

-1 

d 

d 

0 

d 

d 

0 

d 

d 

0 

u 

d 

-1 

d 

u 

1 

d 

d 

0 

d 

u 

1 
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TABLE II. The number of cases N^{Qi 2 , Q 34 ,) and the prob¬ 
ability P^{Qi 2 , Q 34 ) of the charge configurations for two ad¬ 
jacently produced mesons in a flux-tube fragmentation. 


TABLE IV. The number of cases {Q 12 , Q 34 ) and the prob¬ 
ability {Qi 2 , Q 34 ) of all charge configurations for two non- 
adjacent mesons P 42 and P34 in a flux-tube fragmentation. 



Q34=-1 

Q 34 = 0 

T—1 

II 

0 “ 



1—1 

I 

II 

CO 

o> 

Q 34 = 0 

Q34 = 1 


Qi2=-1 

0 

1 

1 



Qi2=-1 

1 

2 

1 

N^{Qi2,Q34) 

Qi2= 0 

1 

2 

1 


N^{Q42,Q34) 

Qi2= 0 

2 

4 

2 


Qi2=+1 

1 

1 

0 



Qi2= + 1 

1 

2 

1 


Qi2=-1 

0 

0.125 

0.125 



Qi2=-1 

0.0625 

0.125 

0.0625 

P^{Qi2,Q34) 

Qi2= 0 

0.125 

0.250 

0.125 


P^{Qi2,Q34) 

Qi2= 0 

0.125 

0.250 

0.125 


Qi2=+1 

0.125 

0.125 

0 



Qi2= + 1 

0.0625 

0.125 

0.0625 


We next examine the production of two mesons in 
the nonadjacent state. We construct in Table |III| all 
charge configurations of two nonadjacent mesons. In 
that case, the charge and flavor quantum numbers of 
quark p 2 and antiquark need not be correlated. With¬ 
out the constraint on the flavors of p 2 and p^, one has 
altogether 16 possible cases as listed in Table III The 
number of cases A^^(Qi 2 ,Q 34 ) and the associated prob¬ 
abilities Q 34 ) in different charge configurations 

(Qi 2 )Q 34 ) for two nonadjacent mesons in a flux tube 
fragmentation are given in Table IV In this case, there is 
an equal probability for two nonadjacent mesons Pi 2 and 
Pm to have charges of the same sign or opposite signs. 
This is in contrast to adjacent mesons for which there is 
no probability that their charges will be of the same sign. 


TABLE III. Quark and antiquark configurations for two non¬ 
adjacent mesons in flux tube fragmentation where P 2 and ps 
are not constrained by flavor balance and local charge conser¬ 
vation. 


Pi 

P2 

Qi2 

P3 

Pi 

Q 34 

u 

u 

0 

u 

d 

-1 

d 

u 

1 

u 

d 

-1 

u 

u 

0 

u 

u 

0 

d 

u 

i 

u 

u 

0 

u 

d 

-1 

d 

d 

0 

d 

d 

0 

d 

d 

0 

u 

d 

-1 

d 

u 

i 

d 

d 

0 

d 

u 

i 

u 

u 

0 

d 

d 

0 

d 

u 

i 

d 

d 

0 

u 

u 

0 

d 

u 

i 

d 

u 

i 

d 

u 

i 

u 

d 

-1 

u 

d 

-1 

d 

d 

0 

u 

d 

-1 

u 

d 

-1 

u 

u 

0 

d 

d 

0 

u 

u 

0 


The charge pattern for two produced mesons in the 
adjacent state is quite different from the charge pattern 
of mesons in the nonadjacent state, in a flux-tube frag¬ 
mentation. 


VII. CHARGE AND FLAVOR CORRELATION 
IN FLUX-TUBE FRAGMENTATION FOR 
QUARKS WITH THREE FLAVORS 


The considerations in the last section for the produc¬ 
tion of u-u and d-d pairs can be generalized to the case 
with the additional production of s-s pairs. Compared 
to the production of a m or d quark pair, the probability 
for the production of a strange-quark pair is suppressed 
by a strangeness suppression factor fs that is energy de¬ 
pendent [46H4^ IMUHT] . and is of order 0.10 at y/s^= 
200 GeV as shown in Appendix A. We can enumerate all 
charge configurations in a flux-tube fragmentation, for 
both adjacent and nonadjacent mesons. The additional 
strangeness degree of freedom brings in an additional cor¬ 
relation of the strange quantum numbers of the produced 
mesons. 

In enumerating all possible configurations for the pro¬ 
duction of two mesons, it is useful to associate each con¬ 
figuration by the order of the strangeness suppression 
factor fs- We use the light q-q pair-production probabil¬ 
ity as a unit of measure of order 1. Upon referring to the 
constituents in Fig. for adjacent mesons and Fig. for 
nonadjacent mesons, the production of an s in pi or an 
s quark in p 4 will each bring in a suppression factor of 
order fs . Similarly, the production of an s in p 2 and an s 
in p 3 in two nonadjacent mesons will also each bring in a 
factor of fs- However, the production of an s in p 2 and an 
s in p 3 in adjacent mesons will bring in altogether only a 
single factor of fs because this s-s pair arises from a single 
Schwinger pair-production mechanism. From these con¬ 
siderations, we can construct all possible meson charges 
(Qi 2 ,Q 34 ) and meson strangeness (5'i2,S'34) up to first 
order in fs in Table V for adjacent meson pairs. Here, 
Sij is the strangeness of the composite meson Py with 
constituents pi and pj. 
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TABLE V. Quark, antiquark, charge, and strangeness config¬ 
urations for two adjacent mesons in a flux-tube fragmentation, 
where p 2 and ps are constrained by charge and flavor balance, 
for flux-tube fragmentation with three flavors. 


Pi 

P2 

Qi2 

S 32 

P3 

Pi 

Q 34 

534 

order 

u 

u 

0 

0 

u 

d 

-1 

0 

1 

d 

u 

1 

0 

u 

d 

-1 

0 

1 

u 

u 

0 

0 

u 

u 

0 

0 

1 

d 

u 

1 

0 

u 

u 

0 

0 

1 

u 

d 

-1 

0 

d 

d 

0 

0 

1 

d 

d 

0 

0 

d 

d 

0 

0 

1 

u 

d 

-1 

0 

d 

u 

1 

0 

1 

d 

d 

0 

0 

d 

u 

1 

0 

1 

u 

s 

-1 

-1 

S 

d 

0 

1 

/. 

d 

s 

0 

-1 

s 

d 

0 

1 

/. 

u 

s 

-1 

-1 

s 

u 

1 

1 

fs 

d 

s 

0 

-1 

s 

u 

1 

1 

fs 

s 

u 

1 

1 

u 

d 

-1 

0 

fs 

s 

u 

1 

1 

u 

u 

0 

0 

fs 

u 

u 

0 

0 

u 

s 

-1 

-1 

fs 

d 

u 

1 

0 

u 

s 

-1 

-1 

fs 

s 

d 

0 

1 

d 

d 

0 

0 

fs 

s 

d 

0 

1 

d 

u 

1 

0 

fs 

u 

d 

-1 

0 

d 

s 

0 

-1 

fs 

d 

d 

0 

0 

d 

s 

0 

-1 

fs 


From the above Table V, we can construct the num¬ 
ber of different charge configurations in Table [VT| for two 
adjacent mesons in a flux-tube fragmentation with three 
flavors. Here, the superscript (ud) specifies the ud sector, 
and (s) the strange sector. 


TABLE VI. Up to order fs, the number of cases 
Q 34 ) and N^^‘\Qi 2 , Q 34 ) for different charge 
configurations (Qi 2 ,Q 34 ,), in two adjacent mesons P 12 and 
P 34 in a flux-tube fragmentation with three flavors. 



Q34=-1 

Q 34 = 0 

t-H 

II 

CO 

0“ 


Qi2=-1 

0 

1 

1 


Qi2= 0 

1 

2 

1 


Qi2= 1 

1 

1 

0 


Qi2=-1 

0 

2fs 

fs 

N^^'^\Qi2,Q34) 

Qi2= 0 

/. 

3fs 

2fs 


Qi2= 1 

2fs 

fs 

0 


In a flux-tube fragmentation, there is a local conserva¬ 
tion of strangeness or heavy-quark flavor in the produc¬ 
tion of a Q-Q pair in the flux tube, similar to the conser¬ 
vation of charge. Hence, there are also correlations that 
arise from strangeness and heavy-quark conservation in 
a flux-tube fragmentation. We present in Table VH the 
number of cases for the production of strangeness con¬ 
figurations (S'12, <S'34), in a flux-tube fragmentation with 
three flavors. 


TABLE VII. The number of cases N^^"‘^\Si 2 , S 34 ) 
and 534 ) for different strangeness configurations 

(5i2, > 534 ), for two adjacent mesons P 12 and P 34 in a flux-tube 
fragmentation with three flavors. 



S'34=-l 

0 

II 

T-1 

II 

Co 

N^^‘’\Si 2 , 834 ) 

Sl2=-l 

834 = 0 
S'l2= 1 

0 

4/. 

0 

0 

0 

4fs 

"^00 


From the above Table lYg the probability 
P'^(“‘^’®)((5 i2, Q34) for the occurrence of adjacent 
non-strange mesons (ud) and strange mesons (s), with 
charges (Qi2,Q34) is related to the number of cases 
jyA{ud,s) Q34) in flux-tube fragmentation with three 
flavors by 


pAM,s)(gi2^Q34) 


^AM.s)(Qi^^Q34) 

8 + 12 /, 


(47) 


which is normalized to 


^ = 1. (48) 

Qi2,Q34 


Similarly, if one considers the strangeness config¬ 
urations of two adjacent mesons, the probability 
P'^(“‘*’®)(S'i2, S'34) for the occurrence of strange mesons 
with strangeness (S'i2,S'34) is related to the number of 
cases iV"^("‘^’'*)(5'i2, <5'34) in flux-tube fragmentation with 
three flavors by 


pA(s)fq q \ (812,834) 

P (^12,834)= ■ 


(49) 


In a similar manner, we can investigate all possi¬ 
ble quark, antiquark, charge, and strangeness configu¬ 
rations for two nonadjacent mesons Pi 2 and P 34 in a 
flux-tube fragmentation. In that case, the charge and 
flavor of P3 will not need to be correlated with those of 
P2- The details of the enumeration is presented in Ap¬ 
pendix B. Again, we need to keep track of the order of 
the strangeness suppression factor /, in different con¬ 
figurations. Keeping configurations only up to the first 
order in /,, we list as a summary the number of cases 
N^^'^d) aji(j 7v^(s) for different combinations of charge 


and strangeness configurations in Tables VHI and [Irre¬ 
spectively. 


TABLE VIII. Up to order /,, the number of cases 
jYV(tid)(gi 2 , N^^'"\Qi 2 , Q 34 ) of the charge configu¬ 

rations for two nonadjacent mesons Pi 2 and P34 in a flux-tube 
fragmentation with three flavors. 



Q34=-1 

0 

II 

CO 

o* 

T-1 

II 

CO 

0“ 


(3i2=-1 

1 

2 

1 

V'^'“‘'HQi2,Q34) 

Qi2= 0 

2 

4 

2 


Qi2= 1 

1 

2 

1 


Qi2=-1 

2fs 

4/= 

2fs 

N’^^‘’\Qi2,Q34) 

Qi2= 0 

4/. 

8fs 

4/. 


Qi2= 1 

2fs 

4/. 

2fs 
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TABLE IX. The number of cases and 

534 ) of the strangeness configurations of two non- 
adjacent mesons in a flux-tube fragmentation with three fla¬ 
vors. 



S34=-l 

534 = 0 

534 = 1 

N^^^\Si2,S34) 

5i2=-1 
Sl2= 0 
5i2= 1 

00 

o 

8/, 

0 

8fs 

00 

O 


From the above Table VIII we can sum up all 
possible numbers of cases, and get the probability 
<534) for nonadjacent mesons. The prob¬ 
abilities for the occurrence of nonstrange mesons (ud) 
and strange mesons (s), with charges (Qi2,Q34) for a 
flux-tube fragmentation with three flavors is related to 
the number of cases (Q12, Q34) by 

Q34) = 


which is normalized to 


^ [P^(“")(Qi2, Q34) + ^"^(*^^12,^34)] = 1 . ( 51 ) 

Qi2,Q34 


We also have the probability {S12, S34) for 

nonadjacent mesons with strangeness ( 5 'i 2 , 5 ' 34 ) for flux- 
tube fragmentation with three flavors given by 




iV^(“^’^)(gl2,g34) 
16 + 32 /, 


( 52 ) 


VIII. TWO-MESON CORRELATION 

FUNCTION IN DIFFERENT CHARGE AND 
STRANGENESS GONFIGURATIONS 

The results in the last few sections allow us to provide 
different signatures for a color flux-tube fragmentation. 
One measures the A(j)-Ar] correlations for a given charge 
or strangeness configuration iz within a pt domain of in¬ 
terest, the pattern of the correlation will reveal the na¬ 
ture of the reaction mechanism within that pT domain, 
as in the case for the hard-scattering process for ^^> 0.5 
GeV/c in Fig. 1 . 

To specify a charge or strangeness configuration v, 
one can choose a pair of mesons with electric charges 
of the same sign, for which (5i2=(534, of opposite sign, 
for which Qi2=~Q34, with the same strangeness, for 
which 812=834, with opposite strangeness, for which 
512=—S'34, or the correlation of a strange meson with 
a nonstrange meson. Different choices will yield different 
probabilities P^ and P^ for finding adjacent or nonadja¬ 
cent pair of mesons, and will lead to different patterns of 
the correlation function, in the fragmentation of a color 
flux tube. 

We can write down the two-meson distribution func¬ 
tion dN/dA(j)dAr] for a pair of mesons with azimuthal 
angle difference Acj) and pseudorapidity difference Arj in 


a charge configuration (or strangeness configuration) v 
within a given px domain in flux-tube fragmentation as 
given by 


dN{u) 

dAcjydAr] 


dN^ A dN^ 
dA(j)dArj dA(j)dArj 


X P^{v). 


( 53 ) 


Here, the first factor in each term, dN^/dA(j)dAr], is 
the distribution function as a function of Acj) and Ary, 
when the meson pair is an adjacent pair with X=A or 
a nonadjacent pair with X=N, in the fragmentation of 
the flux tube, as given by Eqs. ( 44 ) and ( 46 ). The second 
factor P^{v) in each term specifies the probability for the 
pair of mesons in the X state to be in the configuration 
V of interest. The total probability P^{v) for the charge 
configuration v in meson pairs in the state X is given by 


(Qi2 ,< 534 )^ 1 ^ 

( 54 ) 

The total probability P^{v) for the strangeness configu¬ 
ration V to be in the state X is given by 


E [P'"'"'^KS42,834) + P^^^\842 , 834)]. 

(Si 2 ,>S' 34 )£i^ 

( 55 ) 


The quantities ^34) and Q34) 

are given by Eqs. ([ 47 |-([^ and Tables 

We can consider the configuration v of two mesons with 
charges of opposite signs. We obtain from Table |V^ for 
adjacent mesons and Eq. (IrT 


P ^(opposite charge) = =^, ( 56 ) 

and from Table |VIII| for nonadjacent mesons and Eq. 

(opposite charge) = ^ = ^. ( 57 ) 


Similarly, for the charge configuration d of two mesons 
with charges of the same sign, we obtain from Table |VI| 
for adjacent mesons and Eq. ( 47 ) 


P"^(same charge) = 0 , 


( 58 ) 


and from Table |VIII| for nonadjacent mesons and Eq. 

P~(8.me charge) = = 1. (59) 


These probabilities for different charge configurations are 
independent of /,, indicating that they can be obtained 
by considering either two flavors or three flavors. On the 
other hand, the probability for different strangeness con¬ 
figurations must be obtained by considering three flavors. 

We list below the probabilities {P^{u), P^{v)) for a 
few two-meson configurations v. 
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1 . i/=(two mesons with charges of opposite signs), 

(P^, P'^)(opposite charge) = (0.25, 0.125). (60) 

2 . i/=(two mesons with charges of the same sign), 

(P"^, P^)(same charge) = (0, 0.125). (61) 

3 . i/=(two charged mesons {Q 12 ^ 0 and (334 ^ 0 )), 

(pA,pAr)(aii charge) = (0.25, 0.25), (62) 

which is the sum of the first two cases. 

4. :/=(two mesons with opposite strangeness). 


(P"^, P^)(opposite strangeness) = (- 


4 /. 


,0). (63) 


'8 + 12 /« 

5. z/=(two mesons with same strangeness), 

(P'^,P^)(same strangeness) = (0,0), (64) 

which is true only to order /g. The probabilities 
will be non-vanishing to the next order of /^. 

6 . i/=(one |*S '|=1 meson and one nonstrange meson), 

(P"^,P^)(one strange, one nonstrange) 

(65) 


. 8 /g 32 /g 

I lo-f ^ a I 00 ^ 


^8 + 12/g’ 16 + 32/g 


IX. NUMERICAL EXAMPLES OF 
THEORETICAL dN/dA(j)dAri IN FLUX-TUBE 
FRAGMENTATION 


Using Eqs. (44), (46), and (53), we obtain the correla¬ 
tion function for two mesons in a charge or strangeness 
configuration v as given by 


dN^ 


dA(j)dAri M-|- 


1 


+P^{iy)- -^ 


1 


A?)|)/a §7^2 ■ 


( 66 ) 


where dN^/dAcj) is given by Eq. (38) and its numerical 
values are shown in Fig. 

In typical experimental measurements of the two- 
hadron angular correlations such as those in the 
NA61/SHINE Collaboration one collects the 

data for a set of events with similar characteristics (such 
as charge multiplicities) and picks a pair of mesons of 
configuration v (for example, of opposite charges) from 
the same event to build up the “sibling” distribution 
dNsih{v)/dA4)dAri per meson pair, as a function of Acf) 
and A 77 . To eliminate systematic errors and to account 
for the phase space boundaries of the detectors, one also 


uses this set of data to pick a pair of mesons of configura¬ 
tion V from two different events (mixed events) to build 
up the “mixed” distribution dN^iy^{v)/dAtjydAri per me¬ 
son pair as a function of A<j) and A 77 . The sibling distri¬ 
bution contains all the correlation information while the 
meson pair from different “mixed” events do not contain 
any correlation. We can construct the sibling-to-mixed 
ratio 


dNsih{v) / dAtpdArj 
/dA(j)dAr] ’ 


(67) 


which can be compared with theoretical predictions 
C{A(j), Arj) for such a ratio. 

The numerator of (67) can be identified with the dis¬ 


tribution of Eq. ( 66 ). The denominator from different 


events in an event mixing, can be considered in a hypo¬ 
thetical case if the meson pair of configuration v arises 
in the absence of all correlations. With the absence of 
correlations in cfi, then from Eq. (23) dN^/dA(j)=l/8TT'^, 
and with the absence of space-time-ordering correla¬ 
tions in ? 7 , then dNjj/dAr]=dN^/dAr]+dNj^/dAr]=l. 
The distribution dN^iy^{v)/dA(f)dArj in the absence of 
all correlations is [P\v)-\-P^{v)]/8'n'^. It is there¬ 
fore useful to divide the distribution ( | 66 | by the scale 
[P^v)+P^(i/)]/8tt‘^ to construct the following theoreti¬ 
cal angular correlation function for comparison with ex¬ 
perimental ratio of Eq. (67) 


C 


single flux tube 


{A(j),Ar]) = 


dNjiy) 

dA(j)dAri 


[P»+P^(z2)]g^’ 


( 68 ) 


where dN(v)/dA(j)dAr] is given by Eq. ( 66 ). 

It is however necessary to make an amendment in the 
above equation in order to compare with the flux-tube 
fragmnentaion in pp collisions. The above result gives 
theoretical correlation function for a single flux tube. In 
flux-tube fragmentation in pp collisions, two flux tubes 
are formed by the diquark of one nucleon forming a flux 
tube with the valence quark of the other nucleon, and 
vice versa. The mesons from one flux-tube fragmenta¬ 
tion will not be correlated with the mesons from the other 
flux-tube fragmentation, but they are included in the sib¬ 
ling and mixed distributions in the experimental mea¬ 
surement. These uncorrelated mesons from two different 
flux tubes contribute P^{v)/8tt'^ to both the sibling and 
the mixed distributions. Upon taking into account this 
contribution from two different flux tubes, the proper 
theoretical angular correlation function for the charge or 
strangeness configuration v in pp collisions is 




dNi 




C{Acj>,Ap) = 


l-|-e(|AT)|-m)/a dAtj> ' l-|-e(™-|A.7|)/a 


.P^M 


[p» + p^(^)]g^-fp^(^)g^ 


(69) 


which is to be compared with the experimental sibling- 
to-mixed ratio of Eq. (67). 

We can consider an explicit numerical example. For 
pp collisions at ^s,y„= 200 GeV, we have dNch/df]\rj^o ^ 
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2.25 [Bm 66 | . As pp collisions contains contributions 
from two quark-diquark strings, and each string produces 
dNch/dri^l.l2^, we therefore have w~0.59 at y/s = 200 
GeV. 

We can estimate the diffuseness parameter a from the 
ratio of the width (represented by the standard devia¬ 
tion Ut) of the distributions of the pair production proper 
time to the average pair production proper time {r-p^o)- 
Previously in the scatter plot of the pair production 
proper times in the Lund model (see Fig. 2 of m), 
we find 


(7t- 1 

(Tpro) 2 ■ 


(70) 


From Eqs. (39) and (41), w is inversely proportional to 
the pair production proper time Tpro- Hence, we have 


(width of w) (width of Tpro) 

w (Tpro) 


(71) 


The width of w can be identified with 1.5a, for which the 


Fermi distribution of Eq. (43) drops from 82% to 0.18% 


at w—1.5a to ra-|-1.5a. We obtain then 


1.5a Ur 

W (Tpro) 


1 

2 ’ 


and a ^ 0 . 2 . 


(72) 


For numerical purposes, we shall use w = 0.59 and a=0.2 
in subsequent calculations. 


A. Mesons with Opposite Charges 


For the case of a pair of mesons with opposite 
charges, the probabilities for the pair to be in the ad¬ 
jacent and nonadjacent states, as given by Eq. (60), are 
(P^,P^)=(0.25,0.125). 



FIG. 6. (Color online) The two-meson angular correlation 
C(A</), Ari) for two oppositely charged mesons in a flux-tube 
fragmentation in pp collisions. 


In Fig. the contributions in Eq. (69) combine 
together to yield the two-meson angular correlation 
C'(A^, A 77 ), in a flux-tube fragmentation for two mesons 
with opposite charges, calculated for the sample case of 
161=0.59 and a=0.2 in pp collisions. We observe that the 
correlation is suppressed at (A^, Ary) ^ 0 but is signif¬ 
icantly enhanced at {A(f)^Tr, Aij^O). Such a correlation 
arises from the conservation of momentum, charge, and 
flavor in the production of a quark and an antiquark that 
become constituents of the two adjacent mesons. For two 
adjacent mesons that arise from the nearest neighbors, 
their momenta are likely to be azimuthally correlated 
to be back-to-back because a constituent of one of the 
mesons and a constituent of the other meson share the 
same pair-production process at the same spatial point 
that requires local conservation of momentum. Their 
charges are likely to have opposite signs because the con¬ 
stituent of one of the mesons and the constituent of the 
other meson share the same pair-production process that 
requires local charge conservation. 


B. Mesons with Charges of the Same Sign 


For the measurement on a pair of mesons with charges 
of the same sign, the probabilities for the pair to be in the 
adjacent and nonadjacent states, as given by Eq. (61), 
are (P'^,P^)=(0,0.125). Because = 0, the correla¬ 
tion is independent of A(j). 


meson pair with same charge 



FIG. 7. (Golor online) The two-meson angular correlation 
C{A4>, Arf) for meson pairs with charges of the same sign in 
a flux-tube fragmentation in pp collisions. 


We show in Fig. the correlation function C{Acj), Arj) 
in a flux-tube fragmentation for a pair of mesons with 
charges of the same sign. We observe a suppression at 
Ar] ~ 0 indicating that few, if any, pairs of mesons can 
be produced with charges of the same sign at A?; ~ 0 . 
Such a suppression arises because Arj ^ 0 signals the 
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meson to be adjacently produced but mesons produced 
adjacently are forbidden to have charges of the same sign, 
as indicated in Table ED 


C. Charged Meson Pair 


In another example of charge combinations, one de¬ 
tects a pair of charged mesons and measures the 
correlation between one charged meson with another 
charged meson. Then from Eq. (62), we have 
(P"^, P^)=(0.25,0.25). There is an equal probability for 
the two mesons to be an adjacent pair or a nonadjacent 
pair. 


charged meson pair 



FIG. 8. (Color online) The two-meson angular correlation 
C{A4), Ar]) for charged meson pairs in a flux-tube fragmenta¬ 
tion in pp collisions. 

We show the two-meson angular correlation 
C{A(j),Ar]) for two charged mesons for the case of 
flux-tube fragmentation in Fig. It is in fact related to 
the sum of the correlation functions of Figs.[^and[^ The 
suppression at {A(f>, Ar])^0 arises from the suppression 
of charges both of the same sign and opposite signs at 
{A(j), Ar])^0, while the enhancement at (A(()~7r, Ap^O) 
arises from the back-to-back enhanced correlation of 
charges of opposite signs in adjacent mesons. 


D. Mesons with Opposite Strangeness 


F or a pair of mesons with opposite strangeness, Eq. 
(63) gives the probability P^ = /s/(2 -|- 3/s) for the 


pair to be in the adjacent state. The probability 
is zero, up to the first order in /g. It becomes nonzero 
only when we include contributions up to the next higher- 
order terms in //. 

As an example, we examine the case for pp collisions 
at ^SNj\r=200 GeV for which /s=0.1, as given in Ap¬ 


pendix A. Equation (63) then gives (P"^, P^)=(0.043,0). 
We show in Fig. I^the correlation function for two mesons 
with opposite strangeness quantum numbers. The corre¬ 
lation pattern of such a case is quite distinct. It is sup¬ 
pressed at {A(j), Ar])^ 0 and enhanced at {Acj) ^ tt, Arj^ 
0). As one notes from Eq. (69), because is zero, the 
distribution C(A(/, Arj) for this case is indepedent of /g. 


meson pair with opposite strangeness 



FIG. 9. (Golor online) The two-meson angular correlation 
C{A4>, Arf) for two mesons with opposite strangeness in a flux- 
tube fragmentation in pp collisions. 


E. Production of Mesons with Other Strangeness 
Configurations 


For the case of two mesons with the same nonzero 
strangeness, the probabilities for adjacent and nonadja¬ 
cent mesons is zero up to order /g. They are nonzero 
only when we include additional contributions from the 
higher order of order //. 


From Eq. (65), the case of the correlation of one meson 
with strangeness |S'|=I and another nonstrange meson 
with S = Q gives (P^, P^)=(0.086,0.167). We show the 
correlation function in Fig. |10| which has a shape that is 
different from other configurations. 
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one ISI=1 and one non-strange meson pair 



FIG. 10. (Color online) The two-meson angular correlation 
Ari) for one meson with strangeness 151=1 and another 
nonstrange meson with 5 = 0 in a flux-tube fragmentation in 
pp collisions. 


X. COMPARISON WITH EXPERIMENT 

In a pp high-energy collision, most of the produced 
mesons are pions, with some fractions of strange parti¬ 
cles and higher resonances. The strangeness fraction is 
indicated by the Ar+/7r+ ratio that depends on the colli¬ 
sion energy. It is of order 0.1 at y/sNN= 200 GeV and it 
decreases as the collision energy decreases IB1H57] . 

The resonance fraction, relative to the number of pions, 
is of order exp{ —(mresonance —tUTr)/?).®}, and it increases 
with increasing collision energy. For an effective temper¬ 
ature of order 200 to 300 MeV for pp collisions at 
6—200 GeV, the resonance fraction for p or uj mesons is 
of order e~^ to e“^, or 5 to 14 percent. A resonance 
fraction of such an amount would modify and 
by about 10 to 15% and would not change greatly the 
gross features of the two-hadron correlation patterns in 
flux-tube fragmentation shown in the last section. 


A. Comparison with STAR Two-Hadron 
Correlation Data at ^snn = 200 GeV 


The Star Collaboration measured the angular cor¬ 
relation of produced charged hadrons represented by 
Aply/p^ef which is related to the correlation function 
C{ A(j),Ar]) of Eq. (69) by 


Ap Ap 

, - — - 

Y Pref Pref 


l} = Po{G(A</>,Ar;)-l}, (73) 


where Ap=p—pref, p is the correlated distribution from 
sibling events, pref is the reference distribution from 
mixed events, and po ~ -^Pref is the two-dimensional 
angular density averaged over the angular acceptance 


{A(j), Ap) prVHS] . Thus, C{A(j)), Ap) has the same shape 
as Ap/^pref and differs by an overall constant factor and 
an offset. 

The Star Collaboration found that if one separates the 
transverse momentum regions by the boundary pTb=0.5 
GeV/c, the correlation Apjy^p)^ for two oppositely 
charged mesons in the domain below pT<PTb, as shown 
in Fig. IJa), is distinctly different from the correlation 
pattern in the domain above PT>PTb, as shown in Fig. 
Bb)- We mentioned earlier in the Introduction that the 
two-hadron correlation pattern of Fig.|l](b) is a signature 
of the hard-scattering process, with the occurrence of a 
jet (minijet) at (A0, Ap) ^ 0, and another back-to-back 
jet at the away side with a ridge at A(/)~ tt. 

In Fig. |l](a) the small narrow peak at {A(f>, Ap)~0 
arises from e“''-e“ pairs produced by photon conversion 
in the detector and is an experimental artifact. After 
removing the narrow sharp peak at (Acj), Ap)^0 from 
our consideration, we find that the theoretical correlation 
pattern for two oppositely charged hadrons shown in Fig. 
1 ^ for a flux-tube fragmentation has the same gross fea¬ 
tures as the experimental correlation pattern obtained by 
the STAR Collaboration in Fig. |T](a), for two oppositely 
charged hadrons with pt<0.5 GeV/c in pp collisions at 
y^SNN=200 GeV. There is a depression at (Ap, A(j))^0 
but an enhancement at (Ap~0,A(/)~7r). The similar¬ 
ity in the gross features in Figs. Ba) and 1^ indicates 
that in the domain with pt<0.5 GeV/c, the dominant 
particle-production mechanism is qualitatively consistent 
with the flux-tube fragmentation mechanism. Our theo¬ 
retical result from a microscopic approach supports the 
earlier suggestion by the STAR Collaboration [1TH15] to 
use the two-hadron correlation of opposite charges as the 
signature for the soft particle-production process, and we 
identify this soft process as flux-tube fragmentation. 

There are other signatures of the flux-tube fragmen¬ 
tation for mesons with charges of the same signs, or for 
all charged meson pairs. The comparison is complicated 
by the presence of Bose-Einstein correlations for identi¬ 
cal bosons. Such a two-body Bose-Einstein correlation 
is an enhanced correlation at {A(j), Ap)~ 0 and it masks 
the large suppression of the correlation arising from flux- 
tube fragmentation. Eor the correlation of hadron pairs 
with charges of the same sign or for all charged meson 
pairs, it is however necessary to subtract the contribution 
from Bose-Einstein correlations of identical bosons before 
one can make a meaningful comparison. Eor direct com¬ 
parisons, correlation of hadrons with opposite charges or 
opposite strangeness may provide the best signatures for 
flux-tube fragmentation. 


B. Comparison with NA61/SHINE Two-Hadron 
Correlation Data at yAW = 6—17 GeV 

The NA61/SHINE Gollaboration has reported the 
angular correlation C{Ap, A(f)) cx dN/dA<j)dAp for two 
hadrons with opposite charges and pt< 1-5 GeV/c, for 
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FIG. 11. (Color online) NA61/SHINE two-hadron corre¬ 
lation data C{Ari, A(j>) |46H49 | for two hadrons of opposite 
charges with pt<1-5 GeV/c, for pp collisions at different en¬ 
ergies (from Fig. 3 of Ref. |46]). From left to right and from 
top to bottom, the collision energies are y/s^^= 6.3, 7.6, 8.7, 
12.3, 17.2 GeV. 


pp collisions at Ar=6—17 GeV shown in Fig. |^B1 - 
HU] . One can make the following remarks: 


1. The experimental correlation patterns for two 
hadrons with opposite charges in Fig. 11 have the 
same gross features as the signature for the flux- 
tube fragmentation shown in Fig.[^ indicating that 
the dominant mechanism of particle production for 
Pt hadrons in the domain with pt<1.5 GeV/c 
is qualitatively consistent with the flux-tube frag¬ 
mentation. The pattern of the NA61/SHINE two- 
hadron correlation was also shown to be qualita¬ 
tively consistent with the EPOS model IHH] . 
Although there are many different processes and 
diagrams in the EPOS model, it is likely that the 
dominant process responsible for such a two-hadron 
correlation pattern is the flux-tube fragmentation 
part of the EPOS model. 


For yG]vw=200 GeV, the transverse momentum 
boundary has been located at pt&= 0.5 GeV/c 
[35]. This shows that as the collision energy in¬ 
creases, the fractional contributions from the hard- 
scattering process increases, and the boundary pxb 
has been found to move from ppb greater than 1.5 
GeV/c for yG)^=6—17 GeV dH to the lower value 
of pt 6=0.5 GeV/c for yw=200 GeV [H]. The 
search for the boundary ppb in an energy scan will 
map out PTbiy/sNN) more precisely as a decreasing 
function of the pp collision energy 

4. There are additional predicted correlations of two 
hadrons with opposite strangeness and the corre¬ 
lation of a strange meson with a nonstrange me¬ 
son, as shown in Figs, and where the two- 
hadron correlation patterns for flux-tube fragmen¬ 
tation are quite distinct and may be tested with 
the NA61/SHINE experimental measurements. 

5. In pp collisions at NA6I/SHINE energies, baryon 

resonances can be produced by exciting the incident 
projectile and target protons. The decays of these 
baryon resonances will lead to a greater number of 
7 r+ as compared to tt~ and are likely to occur in the 
projectile and target fragmentation regions with a 
diminishing influence in the central rapidity region. 
At central rapidity, the contribution of tt"*' from 
baryon decay has been estimated to be about 20% 
at GeV and it decreases rapidly as a func¬ 

tion of collision energies m- Correlations of two 
hadrons with charges of opposite signs as shown in 
Fig.[TT] are likely to receive a greater contribution 
from the fragmentation of a flux tube than from 
the decay of baryon resonances. 


XI. SUMMARY AND CONCLUSIONS 


2. According to Eq. (411, the width in the Ary direc¬ 
tion is inversely proportional to dN/dri\rir~^o. Thus, 
as the energy decreases, dN/dr]\rf^o decreases, and 
one expects the width in Ary to increase. There is 
a hint of such an increase of the Ary width in the 
data of Fig. |TT]but a more accurate determination 
of the Ary width will be needed. 


3. The correlation pattern of Fig. 11 


remains un¬ 
changed by shifting the pr domain up to pt= 1-5 
GeV/c [30]. This means that in the collision en¬ 
ergy range of ^/sNN= 6 to 17 GeV, the domain 
boundary ppb between the flux-tube fragmentation 
process and the hard-scattering process is greater 
than 1.5 GeV/c. The boundary ppi, can be located 
by searching for a change of the two-hadron corre¬ 
lation pattern for oppositely charged hadrons from 
that of the flux tube fragmentation process of Fig. 



The fragmentation of a flux tube is initiated by the 
production of quark-antiquark pairs along the tube. The 
production of the q-q pairs occurs locally. By the prop¬ 
erty of local conservation laws, a quark and an antiquark 
produced at the same point must balance momentum, 
charge, and flavor. Subsequent to its production, a mem¬ 
ber of the produced q-q pair interacts with its neighbor¬ 
ing antiparticles to form a detected meson. Local con¬ 
servation laws will lead to correlations of the adjacently 
produced mesons. On the other hand, the production of 
mesons along a flux tube is ordered in space-time and ra¬ 
pidity. Mesons that are produced adjacent to each other 
are also close in rapidity. Hence, the proximity of rapidity 
with a rapidity difference Ay (or approximately a pseu¬ 
dorapidity difference Arj) closer than a rapidity width 
w=l/{dN/dy) can signal their adjacent origin from flux- 
tube fragmentation and the correlation in momentum, 
charge, and flavor. Using the rapidity ordering as a signal 
for producing adjacent mesons, we find that dN/dA(j)dAri 
correlation probability or equivalently C{A(j), Ary) for two 
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opposite charges or strangeness mesons is suppressed on 
the near side at Ary) ~ 0 but the correlation prob¬ 
ability is enhanced on the away side at (Acj) ^ tt, Ary ^ 
0 ). 

In addition to flux-tube fragmentation, there is the 
hard-scattering process in nucleon-nucleon collisions. 
The correlation patterns of hard-scattering is well known, 
with a peak correlation at (At)), Ary) ^ 0 on the near-side 
and a ridge along the Ary direction at A(j} ~ tt on the 
away side. The use of the signatures for flux-tube frag¬ 
mentation and hard scattering allows one to separate out 
the dominance of each process in different pt domains 
separated by the domain boundary pTb- 

In our comparison of the theoretical results with ex¬ 
perimental data, we found that the gross features of the 
signature of flux-tube fragmentation for two oppositely 
charged hadrons are similar to those of the STAR [4ll - 
H5] and NA61/SHINE angular correlation data 

for two hadrons with opposite charges in the \ow-pT re¬ 
gion, indicating that the dominant particle-production 
mechanism for low-px particles is qualitatively consistent 
with flux-tube fragmentation. 

We note that whereas the boundary pxb that separates 
the flux tube fragmentation from the hard-scattering 
processes in pp collisions lies beyond pTb> 1-5 GeV/c 
at the NA61/SHINE collision energies of i/s)v]v=6—17 
GeV, it shifts to pTb=0A GeV/c at the RHIG energy of 
^/snn=200 GeV. Thus, the boundary pxb has been ob¬ 
served to shift to a lower value of pTb as the collision 
energy increases. Future determination of pTb in the en¬ 
ergy scan of NA61/SHINE will map out the boundary 
function PTbiy/SNN) as a precise function of the collision 
energy -yJs^N- In addition to being an intrinsic phys¬ 
ical property of the pp collision process, the boundary 
function PTb(,\/SNN) separating the two processes in pp 
collisions also may have implications on the early evo¬ 
lution and thermalization of particles, the quenching of 
jets, and the formation of the quark-gluon plasma, in 
high-energy nucleus-nucleus collisions. 

Returning to the decomposition of particle-production 
mechanisms into soft and hard components in pp colli¬ 
sions, our theoretical result from a microscopic approach 
of the flux-tube fragmentation supports the earlier phe¬ 
nomenological suggestion by the STAR Gollaboration in 
Refs. [ITHiS] to associate the correlation of two hadrons 
with opposite charges to the dominance of the soft com¬ 
ponent of particle production in the domain with pt<0.5 
GeV/c at ys]Vv=200 GeV. However, our result from 
a microscopic approach does not support the second 
phenomenological suggestion in Refs. [4TH45] of a one¬ 
dimensional Gaussian of Ary independent of At)). For a 
better accuracy, a distribution such as Eq. (44) would 
be a better representation of the two-hadron correlation 
function of the ‘soft’ component for hadrons with oppo¬ 
site charges than a one-dimensional Gaussian. 


The interference of the Bose-Einstein correlation 
makes it difficult to use the correlations of two hadron 
with charges of the same sign or two charged hadrons 


as signatures of the flux-tube fragmentation process. In 
addition to the the Bose-Einstein correlation of identical 
bosons, there is also the production of resonances that 
will complicate the pattern of two-hadron correlation. 
The untangling of these many different correlations pat¬ 
terns as a function of collision energies will provide more 
precise information on the changing role of the flux-tube 
fragmentation and hard scattering as a function of the 
collision energy. 

Our successes in identifying the signature of the frag¬ 
mentation of a flux tube will provide a window to ex¬ 
amine further the dynamics of the many-meson system 
and the associated many-meson correlations produced in 
the flux-tube fragmentation. The space-time-rapidity or¬ 
dering stipulates a very simple many-meson space-time- 
rapidity correlations. One may wish to find out how a 
chain of many mesons correlate with each other in their 
azimuthal angles and rapidities, after their production 
in a flux-tube fragmentation. An interesting experimen¬ 
tal and theoretical question is the case of an exclusive 
measurement on an event-by-event basis if the momenta 
of all mesons in the flux-tube fragmentation have been 
measured, and whether it is possible to piece together 
different produced mesons in that event to come up with 
the configuration of the whole system at the moment of 
its fragmentation. These and many other questions may 
be opened up for examination upon a successful experi¬ 
mental search for the signature of the flux-tube fragmen¬ 
tation. 


Appendix A: Relation between K'^ ratio and /s 

In pp collisions at ^/snn^ 6—17 GeV, the ratios 
, K~ Itt~ , and I'k~ depend on the collision 
energy and are affected by the rescattering and decays of 
baryonic resonances formed by the inelastic scattering of 
the incident protons. As a consequence, the ratio Ar+/7r+ 
is greater than K~ /tt” [46H481167] . In higher energy pp 
collisions at y/sNN= 200 GeV, the effects of rescattering 
and the decays of baryonic resonances are less important 
in the central rapidity region, and K'^ j'^K~ j' k~ 
0.10—0.12, and 7r+/7r“~l [M]. It becomes appropri¬ 
ate to estimate the strangeness suppression factor 
at ^/swN= 200 GeV using the flux-tube fragmentation 
model. In such a model. Table V gives the number of 
cases with meson charges (Qi 2 ,Q 34 ) 

in various charge states in the nonstrange sector, and 
the strange sector, for a pair of adjacent mesons. We can 
treat the pair production at all other vertices as a copy 
of Table V, which can then be used to calculate the ra¬ 
tio of the single-particle production of 7r+ and K'^. We 
find from Table V the number of cases of nonstrange {ud) 
mesons with a positive charge, 

^ iV^“‘')(Qi2,g34)=4 + 2/,. (Al) 

Qi2,Q34 —+ 
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We find from table V the number of cases of strange (us) 
mesons with a positive charge, 


^ fV(*Hgi2,Q34)=4/,. 

Qi2,Q34— 

Therefore, we have 

(number of positive strange meson) 
(number of positive nonstrange 7r+ mesons ) 

4 /, 

4 + 2// 


(A2) 


jus) 

{ud) 

(A3) 


Experimentally, at ^/s >200 GeV, the STAR Collabora¬ 
tion gives [64] 

Ar+ (number of K^) 

7r+ (number of 7r+) 

which leads to the estimate 


= 0.10 to 0.12, (A4) 


/= 0.105 to 0.126. (A5) 

Although Table V may appear complicated, it is intu¬ 
itively reasonable that Ar+/ 7 r+ should be nearly equal to 
fa in flux-tube fragmentaion because given a u quark, 7r+ 
is the result of a d-d pair production, and K'^ is the re¬ 
sult of an s-s pair production with a reduced probability 
fa- Hence, Ar+/ 7 r+ should be nearly equal to fa¬ 
ll should be pointed out that the strangeness suppres¬ 
sion factor we have estimated is significantly smaller than 
the value of /s=0.3 used in the standard Lund Monte 
Carlo programs [iniii 7 |. The determination of fa in the 
flux-tube fragmentation from the Ar+/ 7 r+ ratio is best 
carried out when the flux-tube fragmentation process can 
be isolated without interference from contributions from 
other mechanisms. However, for pp collisions at energies 
much below /s = 200 GeV, the production is influenced 
by contributions from baryon resonance production in 
which an incident proton is excited to a baryon resonance 
which subsequently emit a meson llTlIlHj. The baryon 
resonances are produced mostly at the projectile and tar¬ 
get fragmentation regions but have contributions at cen¬ 
tral rapidity when the collision energies are low as in the 
y/s ~10 GeV region. The determination of fa in these 
lower-energy region will be affected by the parameters 
specifying baryon resonance production processes. Only 
for high enough energies at central rapidity can the flux- 
tube fragmentation processes be isolated with small con¬ 
tributions from baryon resonance contributions. Future 
extraction of fa for the flux-tube fragmentation around 
-/s = 200 GeV withpr < 0.5 GeV/c will be of interest in 
clarifying the magnitude of the strangeness suppression 
factor fa- 


Appendix B: Charge and Strangeness 
Configurations for Flux-Tube Fragmentation of Two 
nonadjacent Mesons with Three Flavors 

We examine the quark charge and flavor configurations 
in two nonadjacent mesons P 12 and P 34 in flux-tube frag¬ 
mentation with three flavors as depicted in Fig. The 


constituents p 2 and ps are produced at different vertices 
and are not constrained by charge and flavor conserva¬ 
tion. The cases of possible charge and strangeness config¬ 
urations come in three parts. By including only cases up 
to the first order in fa, we list the charge and strangeness 
configurations in which p 2 is a u quark in Table jX] as Part 
I, P 2 is a d quark in Table |XHI| as Part H, and P 2 is an s 
quark in Table jXVTj as Part HI. The corresponding num¬ 
ber of cases N in the nonstrange sector and 
in the strange sector are given in Tables [Xl] |XH| |XIV1 
pTVlprVlTl and |XVHI| 

In a flux-tube fragmentation with three flavors for non¬ 
adjacent mesons P 12 and P 34 in which p 2 is a u quark, we 
can use Table |X| to construct Table IXII as Part I of the 
number of cases Q 34 ) in the strange and 

nonstrange sectors. 


TABLE X. Part I of quark and antiquark configurations with 
three flavors for two nonadjacent mesons P12 and P34 in a 
flux-tube fragmentation with three flavors in which p2 is a u 
quark. 
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TABLE XI. Part I of the number of cases ^ 

for nonadjacent mesons in different charge states (Qi 2 ,Q 34 ) 
in the nonstrange sector (ud) and the strange sector (s), in 
the fragmentation of a flux tube with three flavors in which 
P 2 is an u quark. 
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From Table [^ we can similarly construct Table XH 
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as Part I of the number of cases {S12, S34) for 

nonadjacent mesons P12 and P34 in different strangeness 
states (iS'12, *534), in the fragmentation of a flux tube with 
three flavors in which p2 is a u quark. 


TABLE XII. Part I of the number of cases 834) 

for nonadjacent mesons P12 and P34 in different strangeness 
states (812,834) in a flux-tube fragmentation with three fla¬ 
vors, in which p2 is a m quark. 
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In a flux-tube fragmentation with three flavors for non¬ 
adjacent mesons P12 and P34 in which p2 is a d quark, the 
possible charge and strangeness configurations are listed 
in Table IXlTll as Part II. 


TABLE XIII. Part II of quark and antiquark configurations 
with three flavors for two nonadjacent mesons P12 and P34 in 
a flux-tube fragmentation with three flavors in which p2 is a 
d quark. 


TABLE XIV. Part II of the number of cases 
^iV{„d,«)(g^2, Q34) for nonadjacent mesons in different 
charge states (Qi2,Q34) in the nonstrange sector (ud) and 
the strange sector (s), in a flux tube fragmentation with 
three flavors in which p2 is a d quark. 
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From Table XIII we can similarly construct Table |XV| 
as Part II of the number of cases (S12, S34) for 

nonadjacent mesons P12 and P34 in different strangeness 
states (5'i2, S'34), in a flux tube fragmentation with three 
flavors in which p2 is a d quark. 


TABLE XV. Part II of the number of cases N^i’’\8i2, 834) 
for nonadjacent mesons P12 and P34 in different strangeness 
states (5 'i 2,5'34), in a flux-tube fragmentation with three fla¬ 
vors in which p2 is a d quark. 
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From Table IXlIIl we can construct Table IXlVl as Part 
II of the number of cases Q34) for nonadja¬ 

cent mesons P12 and P34 to be in different charge states 
(Qi 2)Q34) in the nonstrange and strange sectors, in a 
flux-tube fragmentation with three flavors in which p2 is 
a d quark. 


In a flux-tube fragmentation with three flavors in 
which p2 is an s quark, we list all possible charge and 
strangeness configurations for nonadjacent mesons P12 
and P34 in Table [XVl| as Part III. 

TABLE XVI. Part III of quark and antiquark configurations 
for two nonadjacent mesons P12 and P34 in a flux-tube frag¬ 
mentation with three flavors in which p2 is an s quark. 
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From Table IXVil we can construct Table IXVlII as Part 
III of the number of cases <534) for non¬ 

adjacent mesons P12 and P34 to be in different charge 
states {Qi2,Q34) in the strange and nonstrange sectors, 
in a flux-tube fragmentation with three flavors in which 
P2 is an s quark. 

From the above Table |XVI| we can similarly con¬ 
struct Table IXVIIII as Part III of the number of cases 
IV^(“‘^’®)(5 'i2, S'34) for nonadjacent mesons P12 and P34 in 




































































20 


TABLE XVII. Part III of the number of cases 
Q34) for nonadjacent mesons in different charge 
states (Qi 2 ,Q 34 ) in the nonstrange sector (ltd) and the 
strange sector (s), in a flux tube fragmentation with three 
flavors in which p 2 is an s quark. 
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Upon adding the contributions from all contributions 


from Parts I, II and III, we obtain Tables VIII and IX 
Section VII. 
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